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Abstract — The excitation spectra of two-dimensional strongly 
coupled (liquid-like) Yukawa systems under the influence of a 
magnetic field perpendicular to the plane are found to sustain 
additional high-frequency modes at multiples of the magnetoplas- 
mon. These modes are reminiscent of the well-known Bernstein 
modes but show a number of important differences due to the 
strong coupling of the particles. 

The two-dimensional Yukawa Plasma (2DYP) is being ac- 
tively studied as a model for a diverse selection of physical 
systems, among them the grain- sub system of dusty plas- 
mas [1,2]. The interaction potential for the 2DYP is V{r) = 
(5^/(47r£:o) X exp(— r/A^j/r, where r is the particle distance 
and Q their charge. In contrast to the widely investigated One- 
Component Plasma (OCP), the interaction of the 2DYP parti- 
cles is screened. In the case of dusty plasmas, this screening is 
provided predominantly by the ion component of the plasma. 
The 2DYP is described by two parameters, the ratio of the 
Wigner-Seitz-radius and the Debye length, k = a/Xo, and the 
(OCP-)coupling parameter F = Q'^ /{ATrsoakBT), where ksT 
is the thermal energy of the system. Quasi-two-dimensional 
layers of dust grains are routinely generated in experiments, 
and efforts are currently underway to impose magnetic fields 
on these experimental systems. Furnishing the 2DYP with such 
a strong external magnetic field leads to the modified equations 
of motion [r^ = {xi.ifi)]: 
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where m is the particle mass and B = is the magnetic 
field. To describe the magnetized 2DYP, one additional param- 
eter is necessary, the ratio /3 = uoc/oop between the cyclotron 
frequency uoc = QB/m and the plasma frequency uop of the 
dust grains. The wave spectra of unmagnetized and magne- 
tized 2DYP has been previously calculated by theoretical ap- 
proaches and by computer simulation [2], [3], where, however, 
the emphasis has been on the primary dispersion branches. 
These appear in the vicinity of a combination of the plasma 
frequency uOp and the cyclotron frequency uOc- We have solved 
Eq. (1) by molecular dynamics simulations with integrators 
specifically adapted to the magnetic field for a system with 
periodic boundary conditions. From the microscopic particle 
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trajectories one can obtain the wave spectra in the following 
way: The /c-space current operator 

N 

j{k,t) = Y,^jit)eM^kXj{t)]. (2) 
is separated into longitudinal and transverse currents 

N N 

X{k,t) = Vjx exp[ifexj]; r(fe, t) = Vjy ex.p[ikxj], 

(3) 

from which the fluctuation spectra are computed by a temporal 
Fourier transform J^t ' 
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The viewgraph shows combined MD collective excitation 
spectra L{k,uj) + T(k,uj). At constant n = 2 and F = 200, 
a series of increasingly strongly magnetized systems is shown 
(note that the frequency is scaled by the cyclotron frequency 
cjc). In Fig. 1(a), a moderate magnetization of P = 0.2 is 
depicted which demonstrates the well-known features of mag- 
netized 2DYP: The two distinct transverse and longitudinal 
modes from the unmagnetized case are replaced by two modes, 
the high-frequency magnetoplasmon and the low-frequency 
magnetoshear mode. The magnetoshear mode is an acoustic 
mode, whereas the magnetoplasmon has sl k ^ limit of 
00 = cjc- As we increase the magnetic field [Fig. l(b-c)] 
and shift our focus towards the high-frequency part of the 
spectrum (recall that the scaling factor of the frequency axis 
is ujc which increases linearly with additional dispersion 
branches appear in the mode spectrum. They enter into the 
spectrum as a series of higher harmonics and are, in that 
sense, reminiscent of the classical Bernstein modes. Bernstein 
modes are, however, spaced by uJc- It is evident [cf., i.e.. 
Fig. 1(b)], that the high-frequency modes in the strongly cou- 
pled magnetized 2DYP are not a pure magnetic effect and are, 
consequently, spaced by a combination of uJc and the Einstein 
frequency co'£;(r,/^) [i.e. the average one-particle oscillation 
frequency in a 2DYP] . We have found an analytic formula for 
the mode frequencies [4] 
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which agrees very well the exact spectra, see dotted white 
lines in the viewgraph, over a broad range of magnetic field 
strengths [Fig.l] and all combinations of F and n [4, 5]. By 
considering an additional mode damping, we are also able to 
reproduce the damping of the modes towards small values of 
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ka. Recently the existence of these nonUnear modes could 
be confirmed also for strongly coupled Coulomb systems and 
for systems subject to noise typical for experiments [5]. To 
observe these modes in dusty plasmas it is advantageous to 
use small (radius below 1/im) particles [5]. 
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